Spin glass properties of an Ising antiferromagnet on the Archimedean (3, 12^) lattice 
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We investigate magnetic properties of a two-dimensional periodic structure with Ising spins and 
antiferromagnetic nearest neighbor interaction. The structure is topologically equivalent to the 
Archimedean (3, 12^) lattice. The ground state energy is degenerate. In some ground states, the 
spin structure is translationally invariant, with the same configuration in each unit cell. Numerical 
results are reported on specific heat and static magnetic susceptibility against temperature. Both 
quantities show maxima at temperature T > 0. They reveal some sensitivity on the initial state in 
temperatures where the Edwards-Anderson order parameter is positive. For zero temperature and 
low frequency of the applied field, the magnetic losses are negligible. However, the magnetization 
curve displays some erratic behavior due to the metastable states. 
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I. INTRODUCTION 

The problem of spin glass, with quenched disorder and 
frustration as main ingredients, is known to be a chal- 
lenge in statistical mechanics^ In order to reduce its 
complexity, it makes sense to discuss these ingredients 
separately. The aim of this work is to report numerical re- 
sults on magnetic properties of a structure where frustra- 
tion is present for purely antiferromagnetic interaction, 
with only one value of the exchange interaction. The in- 
teraction is limited to nearest neighbors. The structure 
studied is presented in Fig. Q 

Our motivation to discuss this structure is as follows. 
First, the frustration is to be present, and the simplest 
way to achieve it is to put triangles into the lattice. Sec- 
ond, the number of nearest neighbors is to be odd, what 
ensures that for any magnetic state, the energy barrier 
to flip a spin is finite. Then we may expect that the mag- 
netic transition temperature is greater than zero. Third, 
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FIG. 2: A unit cell. 



FIG. 1: The stretched Archimedean (3, 12^) lattice, with the 
triangular lattice as a background. 



we are interested in the ground state degeneracy. We 
look for a lattice where, if external field is zero, a si- 
multaneous flip of some spins does not change the to- 
tal energy. In fact, if the magnetization of the flipped 
group of spins is zero, the flipping does not change the 
energy even in the presence of the applied held. Fourth, 
all distances between the nearest neighbors should be 
equal. Additionally, two-dimensional structures are pre- 
ferred for their simplicity. The structure presented in Fig. 
^fulflUs all these conditions. The last condition is guar- 
anteed by spanning our structure on the triangular lat- 
tice. Actually, this structure is topologically equivalent 
to the Archimedean lattice (3, 12^). Patterns of all eleven 
Archimedean lattices can be found in Ref. !3. Topological 
equivalence means that one lattice can be stretched into 
the other. Magnetism of Archimedean lattices has been 
investigated for some years«^ However, these works deal 
with the isotropic Heisenberg interaction. We are not 
aware about temperature-dependent simulations on the 
(3, 12^) lattice with Ising antiferromagnetism. 
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TABLE I: List of homogeneous ground states, with spins la- 
beled as in Fig. |21 In last column, paths are indicated which 
lead to other states by flipping pairs of spins in all unit cells. 



Ground 








Pairs 


state 


Sl, S2, S3 


S4 




which can be flipped 


A 


-,+,+ 




-,+ 


(s3, S4) to C or (s2, S5) to B 


B 


-,-,+ 




+ ,+ 


(s2, S5) to A or (si, se) to F 


C 


-,+,- 


+ , 


-,+ 


(sijSe) to E or (53,34) to A 


D 


+ ,-,- 


+ , 


+ , - 


(s3,S4) to F or (82,85) to E 


E 


+ ,+,- 


+ , 




{82,85) to D or (si,S6) to C 


F 


+ ,-,+ 




+ , - 


(si, se) to B or (S3, S4) to D 



II. THE GROUND STATE 

The energy of the classical Ising state^ in the presence 
of external magnetic field H is 

e = ^IjY.s,s,~hY,s. (1) 

where 5^ = ±1, J < is the antiferromagnetic exchange 
constant and the first summation goes over all nearest 
neighbor pairs («, j). The degeneracy of the ground state 
energy (GSE) of the structure presented in Fig. ^ termed 
as stretched Archimedean for brevity from now on, can 
be demonstrated easily when we look at the unit cell in 
Fig. 121 There are nine bonds in the unit cell. The ground 
state energy cannot be smaller than it is allowed by the 
frustration, which is unavoidable within the two trian- 
gles. The contribution from the triangles to GSE is then 
2 J per unit cell. Provided that the spins connected with 
other bonds are oriented in antiparallel, we get 5J per 
unit cell for field equal to zero. This is well possible, and 
some the ground states are periodic. One of such states is 
(si, S2, S3, S4, S5, Se) = (-,+,+,-,-,+)• Another peri- 
odic ground state can be obtained by a simultaneous flip 
of spins S3,S4 in each unit cell. More general, we have 
six periodic ground states (si, S2, S3, —S3, — S2, — si); ad- 
ditional condition is that si-|-S2-fs3 = ±1. These states 
are collected in Tab. ^ 

Obviously, many other non-periodic ground states of 
the whole lattice can be obtained, for example from 
(— , -|-, -f , — , — , -|-) if the flip of two spins is performed 
only in a selected number of unit cells. This means that 
the ground state degeneracy increases with the lattice 
size at least as 6 ■ 2^/^, where N is the number of sites. 



III. NUMERICAL RESULTS 

We apply the heat bath Monte Carlo approach^ to find 
the magnetic contribution to the specific heat at zero 
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FIG. 3: Temperature dependence of the specific heat C. We 
show two curves: for saturated (solid line) and random (dot- 
ted line) initial states. At low temperatures, i.e. below 
T = 0.7 [| J|/fefl], the results depend on the initial state. How- 
ever, the main maximum of C{T) is the same for both initial 
states. 



field, i.e. 

C = P{{E^) ~ {E)^), (2) 

where (• • • ) is an average over thermodynamic ensemble, 
and (3 — l/{kBT). In the numerical calculations, the 
thermal average (• ■ • ) is substituted by the time average. 
A lattice of 6 ■ 10* spins is used, with periodic boundary 
conditions. The initial state is of full saturation, i.e. all 
spins equal to -f 1. Alternatively spins are randomly -1-1 
or —1, with equal weights. After about 100 time steps, 
the total magnetization is close to zero, and the system 
is reasonably close to thermal equilibrium, at least for 
high temperatures. We define a time step as one update 
of the whole network. Then, the time average of E and 
E^ is found. The results are averaged over N^-an = 100 
trials. For T less than 0.7 [| Jj/fes] the statistics is better: 
A'run — 10"^. The obtained plot is shown in Fig. 13 

The same algorithm is used to calculate the static sus- 
ceptibility for zero field, i.e. 

X^PH{M')-{M)% (3) 

where M = X^i*^*- Here again, two sets of data x(T) 
are obtained for two different kinds of the initial condi- 
tions, one saturated and one random. The results are 
shown in Fig. ^ As we see, two curves x(^) become 
different below certain temperature. The curve for the 
random initial state is higher. This is so since, for low 
temperature, the algorithm of heat bath leads the sat- 
urated system to one of its ground states, whereas the 
randomness of initial state is to some extent preserved. 

To get more insight into this memory effect, we 
also calculate the thermal dependence of the Edward- 
Anderson order parameter ^ea-"^ In its definition below. 
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FIG. 4: Thermal dependence of the static spin susceptibihty x 
for zero field. The results for the saturated initial state (solid 
line) show that below the maximum, x tends mildly to zero, 
with numerical uncertainties decreasing at low T. The results 
for the random initial state (dotted curve) show irregularities 
below T = 0.7 [|J|/fcB]. In this range of temperature, we 
observe a rapid increase of x below the main maximum, and 
numerical uncertainties increase when T is lowered. For T — 
0.2 [| J|/A:s], the uncertainties are the largest and reach about 
100 percent of the obtained value. 



*EA 



0.8 - 




1 1.5 
T[|J|/kB] 



2.5 



FIG. 5: Thermal dependence of the Edwards-Anderson order 
parameter. The results are practically the same for initial 
random state and initial saturated state. 



the appropriate time average is written explicitly: 
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(4) 



In Fig. Elwe show $£;^(T) obtained from a random and 
saturated states by the heat bath algorithm at zero tem- 
perature. As we see, ^ea starts to differ from zero at 
temperature close to T ss 0.7 [|J|/fcB], where the plots 
on C{T) and xC^) obtained for different initial conditions 
separate. 

For the case of T = 0, we calculate also the spectrum 
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FIG. 6: (a) Energy distribution of metastable states obtained 
for r = from random initial states. The mean energy is 
about 2.5 % higher that GSE=5J per unit cell, (b) Energy 
against the reduced magnetization m = M/N of the same 
metastable states. No visible correlation is found. m{t = 
0) = 0, H = 0, T = 0, = 6 • 10^ Nrun = 10^ 



of the metastable states. These states are obtained from 
a random initial state by means of the heat bath algo- 
rithm in zero temperature limit. The result is shown 
in Fig. infa). It is seen that there is some randomness 
preserved in the metastable states. The mean energy of 
a metastable state is less than three percent above the 
ground state energy. In Fig. ^h), an attempt is pre- 
sented to find a correlation between energy and magne- 
tization of the same metastable states. As we see, there 
is no correlation at all. We note only that metastable 
states with non-zero magnetic moments do exist and can 
produce some contribution to the hysteresis loop. How- 
ever, this contribution is very small, as seen in the next 
figure. 

The same kind of randomness is present in the results 
on the magnetization curve. The curve is shown in Fig. 
[7| For a comparison, we present also the data obtained 
for the square lattice and the triangular lattice. The 
distribution of sizes of spin flips avalanches is shown in 
Fig. |S1 By an avalanche size we mean the number of 
flipped spins at a given fleld. In Fig. |S1 we see three 
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FIG. 7: The magnetization curves for (3, 12^), the square 
lattice and the triangular lattice (solid, dashed and dotted 
lines, respectively) for T = 0. 



maxima, obtained at fields H ^ 0, H ~ 1 [\J\] and H = 
3 [| J|]. The avalanches at if = are just numbers of spins 
which flip when the system passes from random initial 
state to a metastable state, with final energy distribution 
shown in Fig. Efa). These avalanches are the largest, 
as they contain about 23000 flippings. The avalanches 
occurring at field H = 1 [\J\] and H = 3 [\J\] contain 
10450 and 19800 flippings in the average (see Fig. Ela)). 

Above H = 3 [\J\], all spins are saturated. Provided 
that, in the average, the magnetization of a metastable 
state is approximately zero, a half of spins (3 • 10^) are 
to be flipped to saturate the sample. In Fig. |SIb-d) 
we show the same peaks in smaller scales, what enables 
to observe their detailed character. The positions of the 
peaks of the spectrum indicate, that the summarized size 
of avalanches at if = 1 [\J\] and H = 3 [\J\] is 30250 in 
the average. The small amount of difference, here 250 
spins in the average, mean that some spins flip back and 
forth. 



IV. DISCUSSION 

The magnetic properties of the investigated lattice 
shows similarities to the spin glass state. These are: 
maxima of the spin susceptibility x the speciflc heat 
C as dependent on temperature, and the memory ef- 
fect, i.e. the dependence of these observables on an ini- 
tial state below T = 0.7 [\J\/kB]- In this range of T, 
the Edwards- Anderson order parameter $ ea is different 
from zero. We deduce that there, the ensemble average 
cannot be substituted by the time average, and our nu- 
merical results on C and x are not reliable. However, 
the main peaks of C and x between T = 1.0 [| J|/A:b] and 
T — 2.0 [|J|/A:b] do not depend on the initial state. 

On the other hand, the ground state energy is multi- 
ply degenerated. As a consequence, the overlap- of two 



ground states a and 7, defined as 
1 ^ 



(5) 



is different from unity. For example, the overlap between 
states A and B from Tab. IHis qab = 1/3. When the non- 
periodic ground states are taken into account, the overlap 
distribution P{q) is practically a continuous function of 
q. To explain it, let us consider only two periodic ground 
states, A and B. If only those two are possible, an overlap 
between two non periodic ground states a and 7 within 
a unit cell is either 1 or 1/3, for the cells in the same 
or different states, respectively. For any state a, we can 
select out of A^/6 unit cells and form state 7 putting 
them in the same states, as in a. In this case we have 
the overlap probability 



n, + (iV/6 - n,)/3 



(iV/6)! 



ns\{NlQ-ns)V 



(6) 



In fact, there are six periodic ground states, and not 
only two. Information, cells in which ground states can 
be neighbors without raising energy, is given in the last 
column of Tab. ^ If the condition of periodicity is re- 
moved, the number of ground states of a unit cell in- 
creases. In any case, a typical ground state of the whole 
lattice is expected to be not periodic, but disordered. In 
this way, the most simple definition of spin glass^ is true 
for our lattice. In the case of T = 0, the system dynam- 
ics leads to (meta)stable states, which depend on ini- 
tial states. The spectrum of GSE and the magnetization 
curve reveal a random character, which is a consequence 
of the random initial state. 

Numerical results reported above suggest, that the 
temperature of the transition between the paramagnetic 
state and the low-temperature state, which seems to be 
the spin-glass state, is positive. It is an open question, 
how much disorder is needed to reproduce the vanishing 
of the transition temperature, which is the standard re- 
sult for two-dimensional Ising spin glassesi^ With small 
amount of disorder, we expect that the interaction field 
at some sites is zero and these spins fiip freely. Once a 
cluster of these spins spans throughout the lattice, the 
transition at finite temperatures is likely to disappear. 

Concluding, the (3, 12^) Archimedean lattice can be 
useful as a reference example when the question, which 
features of a realistic spin glass are a consequence of frus- 
tration only, is under debate. This structure can be also 
a good starting point for numerical calculations. Last 
but not least, it is much easier to simulate an Ising spin 
glass where the magnetic interaction is homogeneous. 
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FIG. 8; The spectrum of avalanches, (a) Three peaks are obtained for avalanches at three different fields: (b) H = 0, (c) 
H = 1 [\J\], (d) H ^3 [\J\]. This resuh refiects the distribution of the interaction field. T ^ 0, N = 6 ■ lO"*, iV^un = lO". 
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